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Abstract: In this research, we will propose an efficient FDTD (Finite-difference time-domain) 
calculation method in low frequency range such as kHz or MHz frequency range. This method 
utilize transfer function which is modeled by rational function. The transfer function is calculated 
by the ARMA (Autoregressive moving average model) algorithm to reduce calculation time. The 
effectiveness of the proposed method is confirmed by calculating EMF of steady state 1s excited 
in layered dielectric in the vicinity of the WPT (Wireless Power Transfer) system. 
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1. Introduction 
FDTD method 


Electromagnetic field analysis method in the time domain 


Advantage 
Easy modeling , Good PC memory efficiency, 
suitable for vectorization and parallelization 


iL FDTD Application 


EM field analysis of the Products of complex structure and Human body 
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3 
Products of complex structure 
(Smartphone, TabletPC, etc) Human body and the antenna 
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"FDTD Disadvantage 





Aircraft 


Human body WPT system(kHz or MHz band) 


Need a small cell for modeling Low-frequency EM field Application 





3 By the constraint that Cournt's stability condition 


__Ttrequires Long Caleulation time | 
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For example, Calculation of low-frequency EM field problem 
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Cell size : 5mm cubic cell 
Ti : dt=9.5 
ime step : dt-9.5ps Calculated 
-I- By FDTD time[us] 
Low frequency wave , , 
(MHz or kHz band sin-wave) Time until steady state : T=18pus 

Time step : dt=9.5ps 


The number of time steps 


= T/dt 
1,900,000 steps 





Voltage[|V | 





time[us] 
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Purpose of This study 


Electromagnetic field analysis in low frequency 
Analysis of the antenna characteristics 
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Two proposed method 
1. Estimation of Rational functions with the Auto-Regressive Moving Average(ARMA) 





Value of the time domain ===) ` Hy 
x(t) : Input signal Calculated à | or & 
y(£) : Output signal by ARMA ~ $| 
Frequency|Hz| Frequency|Hz] 
Input Impedance Transfer function 


2. Estimation method of Electromagnetic(EM ) field using a transfer function 
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Estimated 
Frequency[Hz] EM field 
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Transfer. functio 





E|V/m] 


ds 












































pulse Transfer function 





Calculation contents in this study 





Feed Point Resistor =50 Q 


WPT system near layered dielectrics 


Tested for two types of problem 


1. Steady state E field in time domain 


2. Reflection coefficient of loop antenna in frequency domain 


Lb We will evaluate 


Calculation accuracy and The number of time steps 1 
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2. Proposed method 
2.1 ARMA/FDTD Algorithm 


To approximate the rational function in the frequency domain 
Y(Z) ag t+ ayZ~* + a2QZ~* ++ tagZ~4 
X(Z) 1-bhZ--4bj2-4--bjZ7" 


p,q : The order of the rational function model 
T : Sampling time 


H(Z) = **(1) Z — exp(joT) 

















"FDTD 
4 Solve the (2) in order to determine 
x 
Equation of ARMA model 5 
p q 5 
yG) 2 — » by(n- D » ai(n-]) **(Q) | 
= = time|s | 
x(n) : Input signal (ampressed voltage or current) -ARMA/FDTD 





y(n) : Output signal (electric field , current, etc.) 


We checked that all poles are placed the stability region in z-plane 
by using unknown coefficients a; and bj. 


output y(t) 











time[s] 
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2.2 Estimation method of Electromagnetic field using a transfer function 























Input signal ARMA algorithm 
> 
T: 3 
S 9 
S 8 
fam, — 
time[ns] = E 
S z, 
Calculated by E zs 
FDID Frequency|Hz| 














Output signal Calculated by time[s] 
ARMA algorithm 
Estimated by 
Transfer function 


E[ V/m] 














time|ns| 





3. Simulation result 
3.1 Calculation of the steady state E field 
Analysis Model : Wireless Power Transfer System (WPT System) 


Calculation of the steady state E field when sin wave of 5.5MHz was the input signal 
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a Resistor =50 Q 
Feed Point 








Calculation point 


WPT System model Skin muscle bone 
skin :¢=547 [-] o=0.122 [S/m] 
muscle : ¢=282 [-] o=0.594 [S/m] 
bone : ¢=55.1 [-] o=0.0380 [S/m] |19 
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"The calculated data for Transfer function 
Calculation of Transfer function 





yen =- LO -FDO y) : EW [V/m] The E field at each calculation point 
Ww X(Z) T[x(t)|] x(n): Input Voltage [V] m | | 
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á e e 520 30 ae 
Comparison Result : Transfer function(skin) aa 
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0.020 i SSS reac Pai ; bone 
Original FDTD-Real = skit muscl 
ARMA/FDTD-Real y Unit : mm 
0.015 Original FDTD-Imag 
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Transfer function Hx 
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-0.010 
-0.015 (D 
0.020 Calculation point 


frequency [MHz] 


Calculation accuracy : Original FDTD = ARMA/FDTD 
The number of time step : Original FDTD » ARMA/FDTD 


(240,000steps) — (80,000steps) A 
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Calculation Result : E field in time domain(x -direction component) 


Original FDTD(skin) 
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Calculation point 
Skin muscle bone 





time [us] 


Original FDTD = ARMA/FDTD 
: Original FDTD > ARMA/FDTD 13 
(1,900,000steps) (80,000steps) 


sss MMI-NI 


Calculated accuracy : 


The number of time steps 





3.2 Calculation of the Reflection coefficient 
Modell : Loop antenna of WPT system 
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Cte 
100 = = 










~ —* 
CMS © © Resistor =50 Q 


Wie Teed Point 


Unit : mm 


Model2 : Loop antenna of WPT system near three layered dielectric 







en Eod Doint Resistor =50 Q 
Unit : mm 


skin muscle bone 


skin : €=547 [-] o=0.122 [S/m| 
muscle : €=282 [-] o=0.594 [S/m] 
bone  :&,755.1 [-] 070.0380 [S/m] 
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"The calculated data for Reflection coefficient 


Calculation of Impedance 
Y(Z) -ly(t)] 


H(Z) — 


X(Z) [x(t)] 
y(n) : Input Voltage [V] 
x(n) : Input Current [A] 


‘Input Voltage Calculated 
output signal y(t) by FDTD 


Gaussian pulse 
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The input current at each model 
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Comparison Result : Reflection coefficient 


Modell : Loop antenna of WPT system 
Model2 : Loop antenna of WPT system near three layered dielectric 


Modell 
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frequency[ MHz] 


Calculation accuracy : Original FDTD = ARMA/FDTD 
The number of time steps : Original FDTD >ARMA/FDTD 16 
(1,000,000steps) (100,000steps) 
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4. Conclusion 


Y Our proposed method 
1. ARMA/FDTD method 
2. The estimation method of EM filed using a transfer function 


Y We tested the our proposed methods 
by the calculation of the steady state E field and Reflection coefficient. 
The simulation result 


Calculation accuracy : Original FDTD = Proposed method 
The number of time steps : Original FDTD >> Proposed method 


Our proposed method can reduce calculation time of FDTD method 
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